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The global mass dependence of the nuclear symmetry energy asym {A) and its two basic 
ingredients due to the mean-level spacing e{A) and effective strength of the isovector 
mean-potential k{A) is studied within the Skyrme-Hartree-Fock model. In particular, 
our study determines the ratio of the surface-to- volume contributions to asym{A) to be 
rg/Y ^ 1.6 and reveals that after removing momentum-dependent effects by rescaling s 
and K with isoscalar and isovector effective masses, respectively, one obtains e* n* . 

1. Introduction 

The most common route in constructing effective microscopic nuclear models starts 
from the nuclear equation of state in infinite nuclear matter and use such parameters 
like isoscalar saturation density poi the volume binding energy ay, the incompress- 
ibility parameter K^o, and the asymmetry energy oiQn as primary constraints for 
these models. All these values are extrapolated, however, from the studies of finite 
nuclei, in particular from the semi-empirical mass formulas: 

(/^ + Ai)+,.., (1) 
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where / = \N — Z\/A while as and aiym are coefficients defining contributions from 
the surface energy and the surface part of the nuclear symmetry energy (nse), 
respectively. 

Concerning the isovcctor terms in finite nuclei which arc subject of this work 
there is at present no consensus concerning the magnitude, A, as well as origin of 
the Wigner, ~7, term. The Weizsacker mass formula (ld) sets A = 0. The ld of 
MoUer et al. ^ admits only a volume-type Wigner term, which is not consistent with 
Eq. (1). For this term it gives A « 0.975 i.e. a value close to the one used in shell- 
model inspired mass formulas which are using ~ T(T + 1) (i.e. A = 1) where T = 
\Tz\ = — Z\/2. Another controversy exists concerning the surface contribution to 
the NSE. The values of the surface-to-volume ratio rs/v = cii^m/a'^^m quoted in the 
literature vary strongly. For example, Danielewicz ^ estimates it to be 2.0 < rs/v < 
2.8, the mass formula of Ref. ^ yields rs/v ~ 1-6 while the hydrodynamical-type 
models that include properly polarization of the isovector density predict rs/v ~ 
2 3. 

The main objective of this work is to study various aspects of the nse emerging 
within the Skyrme-Hartree-Fock (she) model. This work supplements our earlier 
studies on this subject ^■^'^■'^ to which we refer our readers for details. 

2. The Skyrme- force induced nuclear symmetry energy and its 
mass number dependence 

In Ref. ^ we have demonstrated that the SHE symmetry energy behaves rather 
unexpectedly according to the formula: 

Ei'y^^^ = y{A,T,)T' + ^K{A,T,)TiT+l) , (2) 

where £{A,Tz) « e{A) and k{A,Tz) ~ i^iA) are fairly independent on T^, at least 
for Tz > 8. Hence, for further quantitative analysis of the mass dependence of 
the NSE we use the mean values of £{A) and k{A). These averages over at fixed 
A are calculated using the following restricted set of nuclei: > 4 for ^4 = 20; 

> 6 for A = 24; and > 8 for 28 < A < 128. By using a restricted set of nuclei 
we smooth out both e{A) and k{A) curves in order to diminish the influence of 
shell structure. The e{A) and k{A) are interpreted as the mean-level spacing at the 
Fermi energy in iso-symmetric nucleus and the effective strength of the isovector 
mean-potential, respectively. 

The relation (2) holds extremely well within the SHE model except for the SkO 
parameterization ^ for which we observe an enhancement in the linear part of the 
second term in Eq. (2). For the SkO the second term becomes ~ kT{T + A) with 
A > 1 , as illustrated in Fig. 1 . The value of A is however much weaker than Xb,m f ~ 
\ + e/R found recently in relativistic mean-field (rmf) ^ . It should be mentioned 
that the SkO is characterized by an exceptionally strong isovector strength of the 
spin-orbit term, inspired by RME, having opposite sign compared to standard SE 
parametrization. 
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Fig. 1. Values o{ K.{A,Tz) calculated using Sly4 (circles) and SkO (squares). The calculated points 
mark the percentage of calculated nuclei for which k(A,Tz) deviates from the mean- value ii(A) by 
less than three percent versus ■ Two different scenarios are assumed for the linear contribution 
to ~ kT{T + X) term with A = 1 (filled symbols) and A(= Xumf) = 1 + £/S (open symbols). Note 
that while for the SLy4 A = 1, for the SkO it is half way between 1 < A < Xhmf- 
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Fig. 2. The isoscalar effective mass scaled values of e*{A) = ^^£(A) (stars) and R(A) (circles) 
calculated using the Shf method with SLy4 parametrization. Open symbols denote £*{A) and 
ii{A) averaged over all the calculated nuclei. Filled symbols mark smoothed values of £*{A) and 
h{A) calculated using a restricted set of data. Vertical arrows indicate major shell gaps. Note the 
strong influence of shell structure on s*{A) and the smooth behavior of k{A). 



The global mass dependence of the two components of the symmetry energy, 

e and R is shown in Fig. 2. The figure reveals several universal features which 
appear to be independent of the type of the SF parametrization including: (i) strong 
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dependence of e{A) on kinematics (shell effects); (ii) almost no dependence of 
k{A) on kinematics; (iii) clear surface (~ -^jw) dependence reducing the dominant 
volume term -^) in both s{A) and k{A). 

Indeed, the values of e{A) show characteristic kinks close to double- (semi) magic 
A- numbers. These kinks are magnified when all the calculated nuclei are used (no 
smoothing) to compute £{A), but without affecting qualitatively the overall profile 
of the curve. On the other hand, k{A) is almost perfectly smooth with barely visible 
traces of shell structure. It confirms our earlier conclusion ^ that the gross features 
of the Skyrme isovector mean potential can be almost perfectly quantified by a 
smooth curve parametrized by a small number of global parameters. 




8r2/r2[%] 

Fig. 3. Histogram showing the deviation between neutron skin thickness estimated using the simple 
hydrodynamical formula (3) taken for rg/y = 1.6 and the results of microscopic SHF calculations. 

Note, that the deviation is of the order of ±3%. Note also, that the ccntroid of the histogram is 
slightly shifted to the right what indicates that rg/y should be slightly larger rg/y ~ 1.65. 



The SHF models yield rg/v ^ 1-6 in accordance with the ld ratio ^. There 
are several observables which are strongly sensitive to rs/v These include in par- 
ticular the static dipole polarizability (sdp) and neutron skin thickness. Using the 

hydrodynamical model of one can derive simple expressions for both the SDP and 
the neutron skin thickness. It appears that for rg/y ~ 1-6 these simple expres- 
sions yield results that are very consistent both with the data (for SDP) and the 
microscopic SHF calculations (neutron skin thickness), see ^. For example, the hy- 
drodynamical formula for the neutron skin thickness [in the lowest order expansion 



in -j] reads: 



(r2) A 1^ + 3^1/3 ■■■(■ ^-^^ 
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The interesting feature of Eq. (3) is that it does not depend on the bulk nse coeffi- 
cient but only on the ratio rs/v This formula can easily be cross-checked with our 
microscopic SHF calculations and the results are depicted in the form of a histogram 
in Fig. 3. 
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Fig. 4. The correlation between the effective mass scaled volume 6y and Ky (open dots) and the 
surface and Kg expansion coefficients. Note, that except for the SIII and SkM* interactions, 
the expansion coefficients are equal. 



The most striking result of our analysis is the near-equality of £* » k* [where 

K* = R— and ml denotes the standard isovectorial effective massl occurrine; for 
all modern parameterizations, see Fig. 4. Indeed, e* differs from k* only for old 
parameterizations like the SIII and SkM*. This result confirms the rather loose 
claims often appearing in the textbooks that "the kinetic energy [efg] and the 
isovector mean-potential contribute to the agym in a similar wa-if is indeed correct 
but only after disregarding non-local effects. To our knowledge, it has never been 
discussed why this apparently independent quantities should be similar. 

3. Summary 

The global mass dependence of the nse strength asym{A), e{A), and R{A) is studied 
in detail within the SHE theory. Our study enables us to establish the surface-to- 
volume ratio of asym(A), rs/v ~ 1-6- This value is in agreement with the ld of 
Ref. ^ and is consistent with simple hydrodynamical estimates for the SDP and 
neutron skin thickness. Our study also reveals an almost linear correlation between 
Cf and (mj)^^, see Fig. 5, and a striking similarity between e* k, R* ''. The latter 
near equality indicates that the contribution to agym due to the mean-level spacing 
and due to the mean-isovector potential are similar but only after disregarding non- 
local effects. Whether this is a fundamental property of the nuclear mean field is 
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Fig. 5. Correlation between and kinematic isovectorial cfTcctivc mass (mj;)~^ = ^^Cfpo for 
different Skyrme force parameterizations. The and are coupling constants defining C^pf 
and CI pin isovector terms in the Skyrme local energy density functional. 

an open question that requires further studies. 
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